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Abstract. We study the most central regions of AGN jets with an angular 
resolution of tens of micro-arcseconds using VLBI at millimeter wavelengths 
(mm- VLBI) . We present and discuss a new 86 GHz VLBI survey of compact 
radio sources. We show new high dynamic range images of two nearby radio 
galaxies (3C 120 and M87). In M87 the size of the compact VLBI core (the 
jet base) is < 15 Schwarzschild radii. Future mm- VLBI observations at 1 mm 
and shorter wavelengths should lead to images of galactic and extragalactic 
radio sources with a spatial resolution down to a few Schwarzschild radii of the 
central super massive black holes. To achieve this, the participation of large and 
sensitive millimeter and sub-millimeter telescopes in VLBI is essential. 

1. Introduction and Scientific Motivation 

Radio galaxies and quasars produce powerful relativistic plasma jets, acceler- 
ating particles up to TeV energies. Jets are not unique to AGN, but are also 
observed in other objects like in X-ray binaries, micro-quasars, young stellar ob- 
jects and probably also in Gamma-ray bursters. In modern tiieory, tiie launcii- 
ing and acceleration of the relativistic jets is explained via MHD processes, with 
magnetic fields either attached to the rotating accretion disk or directly linked 
to t he ergosphere o f the central and possibly rotating super-massive black hole 
(e.g. iNaravan 20051 ). 



In order to better understand the process of jet formation, a study of the 
inner most regions of compact radio sources is necessary. VLBI at short mil- 
limeter wavelengths (mm- VLBI) provides, via its unrivaled angular resolution 
of only a few ten micro-arcseconds (at A = 1 mm), an almost unique tool to test 
and discriminate between the different theoretical models. In addition to the 
high spatial resolution (corresponding to < 10 Schwarzschild radii for nearby 
sources like Sgr A* or M87), the high observing frequency allows to image emis- 
sion regions, which are self-absorbed and therefore not observable at longer 
wavelengths. 

In the following we present some results from global 3 mm- VLBI performed 
with the Global Millimeter VLBI Array (GMVAQ) at 86 GHz. VLBI observations 



^see: http://www.mpifr-bonn.mpg.de/div/vlbi/globalmm/index.html 
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Author Ant. Sensitivity Number Detections 



Beasley et al. 1997 


3 


- 0.5 


N= 


=45 


16% 


Londsdale et al. 1998 


3-5 


- 0.7 


N= 


=79 


14% 


Rantakyro et al. 1998 


6-9 


~ 0.5 


N= 


=67 


24% 


Lobanov et al. 2000 


3-5 


~ 0.4 


N= 


=28 


93% 


Lee et al. 2008 


12 


- 0.2 


N= 


=127 


95% 



Table 1. VLBI surveys at 86 GHz: Col.l gives the reference, col. 2 the 
number of antennas, col. 3 the typical detection sensitivity {7a) in [Jy], col. 4 
the number of observed sources and col. 5 the percentage of detected objects 
(detection rate). 

at shorter wavelengths (A < 3 mm, v > 100 GHz) are possible but are still 
limited by the number of available mm-telescopes (typ ically < 3 — 4), their 
sensitivity and the limited observing bandwidth (see e.g. iKrichbaum et al. 20071 . 
and references therein). To date, it is therefore not yet possible to produce 
a reliable map of a radio source at these shorter wavelengths. The GMVA 
combines the sensitive and mm-capable telescopes in Europe (100 m Effelsberg, 
6x15 m Plateau de Bure, 30 m Pico Veleta, 20 m Onsala, & 14 m Metsahovi) 
with the VLBA. The GMVA is open to the scientific community and calls for 
observing proposals twice per year (proposal deadlines: Feb. I''*, Oct. 1'^*). It 
enables imaging of compact radio sources at an angular resolution of ^ 40/ias. 
The transatlantic detection sensitivity (7cr) ranges between ~ (80 — 200) mJy on 
the most sensitive baselines (for 512 Mbps). This should be compared with the 
~ (300 — 400) mJy baseline sensitivity of the stand-alone VLBA. In a typical 
12 hr GMVA observing run, an image rms of fTmap ^ (1 ~ 4) mJy is reached. 

2. A New and Comprehensive 86 GHz VLBI Survey 



VLBI surveys at the highest possible frequency and resolution are important for 
answering questions on the compactness of extragalactic radio sources, system- 
atic differences between source classes, and if and how the brightness temper- 
ature varies along the jet and towards the region, where the jet is made. The 
latter is closely related to observational tests of jet launching models. Due to 
the limited sensitivity of early 3 mm- VLBI and the relatively small number of 
telescopes, previous surveys revealed only moderately low (< 20 — 30%) detec- 
tion rates and small numbers of observable sources. Table [T] summarizes the 
existing 3 mm VLBI surveys in chronological order. 

After more sensitive VLBI telescopes became available around year 2000, 
a new and larger than previous VLBI snap-short survey was conducted. The 
observations took place in 3 sessions (October 2001, April & October 2002), with 
the participation of the VLBA (8 antennas). Haystack, Pico Veleta, Plateau de 
Bure, Effelsberg, Onsala and Metsahovi. Mainly due to the high sensitivity of 
the IRAM telescopes, a detection sensitivity of ~ (0.2 — 0.3) Jy could be reached. 
This resulted in an impressive detection rate of 95%, which equals to detection 
of 121 out of 127 observed sou rces, of wh ich 109 objects (86 %) could be imaged 
(see Table ^ (|Lee et al. 20071 : iLee 20081 ). 

After fitting of Gaussian components to the visibility data for each source, 
the flux, position and size (FWHM) of the VLBI components were measured. 
In Figure [1] we plot the resulting distribution of the measured core flux densi- 
ties (top left) and sizes (bottom left) for Quasars, BL-Lacertae Objects, Radio- 
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Figure 1. Distribution of flux density (top left) and angular size (bottom 
left) of the core components for the survey. The shading denotes different 
object classes. Left: Binned, averaged and normalized visibility amplitude 
plotted versus uv-distance. Different symbols identify the different object 
types. Numbers in brackets give the number of objects in each class. Some 
visibility b ins (interval: 200 MA) do not contain data for each source class 
(|Lee 2OO80 . 

Galaxies and the unidentified objects. The right panel of Figure [T] shows the 
mean normalized visibility amplitudes of the optically identified sources. 

3. The inner jet of 3C 120 

In Figure [2] we show a 43 GHz VLB A map and a 86 GHz GMVA map of the radio 
galaxy 3C 120, observed in October 2004. The left panel shows uv-tapered maps 
convolved with a (0.1 x 0.3) mas beam. The right panel shows the un-tapered 
central 0.5 mas region of the jet, convolved with the smaller 86 GHz beam (0.06 
X 0.25 mas). At the redshift of 3C 120 {z = 0.033) an angular scale of 1 mas 
corresponds to a spatial scale of 0.63 pc. Assuming negligible variability on the 
16 day interval between the two observing dates, it is possible to determine the 
spectral indices a43/86GHz along the jet {S oc z/°). For the labeled regions A to 
F in Figure [2] (left), we obtain spectral indices of -0.35 for A (r = Omas), -0.07 
for B (r = 0.1 mas), -0.36 for C (r = 0.25 - 0.4mas), -0.46 for D (r = O.Smas), 
-0.7 for E (r = 1.2 mas) and -0.8 for F (r = 1.8 — 2.2 mas). A general trend 
of spectral steepening with increasing core separation is obvious. We note that 
component B exhibits a much flatter spectrum than the core component A. With 
a size of ~ 0.1 mas (FWHM) and a spectral turnover frequency (43 < Um < 
86) GHz, component B cannot be in energy equipartition, and is dominated by 
the magnetic field energy. Similar arguments for the core region (component 
A, size < 32^as < 7 ■ 10^^ cm) and for the outer jet components (C-F) yield 
no evidence for a substantial deviation from equipartition. The higher angular 
resolution at 86 GHz (Fig. [21 right) also reveals some structural difference to 
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Figure 2. VLBI maps of 3C 120 at 43 GHz (top, VLBA: Oct. 28, 2004) and 
86 GHz (bottom, GMVA: Oct 16, 2004). The maps result from a reanalysis 
of previously published data (jMarscher et al. 20071 ). The uv-tapered maps 
on the left show a faint and partially resolved one sided jet extending to at 
least r 2 mas. On the right, the central 0.5 mas region is shown with higher 
angular resolution (beam: (0.06 x 0.25) mas). Lines and labels guide the eye 
and help to identify corresponding emission regions. Circles in the contour 
maps on the 2 right panels show size and position of Gaussian components 
fitted to the visibility data. At 86 GHz, the core size is < 32/ias. 

the 43 GHz map. Remarkable is the difference in the position angle of the 
inner (r < 0.5 mas) jet components. This could be interpreted as evidence for 
transverse jet stratification. 

4. 3C 454.3 after a major outburst in 2005: 

In May-June 2005, the quasar 3C 454.3 {z = 0.859) showed a major flux density 
outburst, which was observed from hard X-ray to radio bands (e.g. Villata et 
al. 2007, and references therein). The flux variations at 1.3mm and 0.85mm 
are shown in Figure [3] (top, right). During the flare maximum quasi-periodic 
oscillations on a ~ 11 day time scale appear, possibly indicating motion in a 
helical jet with ~ 100 Schwarzschild radii diameter (assuming Mbh = 10^ Mq). 
A Doppler-boosting factor of 5 ~ 5 — 9 would be necessary to bring the observed 
variability brightness temperature of Tb — 10^^ K (at 230 GHz) back to the in- 
verse Compton limit of ~ 10^^ K. The spectral evolution of the flare through 
the radio bands is shown in Figure [3] (bottom, left). It is obvious that the vari- 
ability is most pronounced at short mm- wavelengths, with a spectral turnover 
moving from initially ~ 230 GHz to lower frequencies. Below ~ 10 GHz almost 
no variability is seen. 

So far, VLBI observations at frequencies < 43 GHz did not reveal strong 
evidence for the ejection of a new jet component, which could be related to this 
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Figure 3. Left: 43 GHz VLBA (top, data from Marschcr et al.) and 86 GHz 
GMVA (bottom) map of 3C 454.3 observed in Oct. 6 and Oct. 17, 2005. Both 
maps are restored with a common beam of (0.1 x 0.25) mas. Lines with la- 
bels connect the three corresponding emission regions A, B and C (see text). 
Top right: Flux density variability at 230 and 345 GHz. The insert shows th e 
variability during the time of the peak flux (data: SMA. lGurwell et al. 20071) . 
Bottom right: Evolution of the radio spectrum of 3C 454.3 after the optical 
flare. The spectra cover a time range from Sep. 2004 until Nov. 2006 at time 
intervals of a few weeks. The arrows mark the direction of spectral evolu- 
tion. Each spectrum is obtained from quasi-simultaneous flux measurements, 
combining data from Effelsberg, VLA, IRAM, SMA & Ratan-600. 




Figure 4. 86 GHz VLBI images of M 87 (3C 274) obtained with global 3 mm 
VLBI in April 2002, 2003, and 2004 (from left to right). The contour levels 
are at -0.3, 0.3, 0.6, 1.2, 2.4, 4.8, 9.6, 19.2, 38.4, and 76.8 % of the peak flux, of 
0.52 Jy/beam (2002), 0.73 Jy/beam (2003), 0.35 Jy/beam (2004). The lowest 
contour level is omitted in the image of 2002. The restoring beam sizes are 
0.23 X 0.057mas (2002), 0.31 x 0.062mas (2003), and 0.20 x 0.054mas (2004), 
respectively. 
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flare. In Figure [3] (left panel) we show two VLBI maps obtained about 4 months 
after the flare. The 86 GHz data reveal a clear elongation of the core region A, 
which is not visible in the 43 GHz map (region B). A spectral index calculation 
for the three main emission regions yields for the core A a43/86GHz = —0.4, for 
the region B a = +(1.1. ..2. 6) and for the jet C a = —0.3. The strongly inverted 
spectrum of the region B indicates possible fore-ground absorption, perhaps 
from a cloud or torus. With a radial extend ranging from r^i^ ^ 0.1 mas to 
''max — 0.3 mas (0.1 mas corresponds to 0.75 pc) and a sufficiently high opacity, 
the inner portion of the jet would become visible only at the highest frequencies. 



5. The Jet of M 87 



In Figure U] we show 3 images of the inner jet in M 87 (3C 274) obtained with 
the GMVA at 86 GHz during 2002-2004. The map of 2004 has a higher dynamic 
range than the earlier images due to increased observing bandwidth (512 Mbps). 
To our knowledge, these are the highest angular and spatial resolution maps 
ever made of this source. At a distance of D=16.75 Mpc, the minor axis of the 
GMVA observing beam of ~ 54^as translates into a spatial scale of 1.35-10^^ cm 
(4.4 mpc). From the slightly elliptical observing beam one obtains a formal up- 
per limit for the core size of <(200 x 54) //as, which corresponds to Rs={56 x 
15) Schwarzschild radii (for Mbh = 3 • 10^ Mq). We may identify the size of the 
brightest and most compact VLBI component with the jet diameter at or near its 
origin. In this case it is remarkable that the jet base is so small and bright (Tg > 
2 • 10^*^ K), contrary to expectation s from magnetic-slingshot type acceleration 
models (e.g. lBlandford Sz Pavne 1982 ). Here, magnetic flelds are anchored in the 
rotating accretion disk, expand and form a light-cylinder. Its diameter define s 
the minimum jet width of > 50 Rs (|Camenzind 19901 : iFendt k Memola 200lh . 
The observed small size therefore points more towa rds models in which the je t 
is attached more directly to the rotating black hole ( Blandford Sz Znaiek 19771) . 



Mc Kinnev 2006l i 



gaining its energy via M HD interaction wit h the inner disk (e.g. 
or the Penrose process ( Gariel et al. 20071 . and references therein). 

The images of 2002 and 2004 both show a bifurcated conically opening 
one-sided jet, similar to the jet morphology seen at 43 GHz (Walker et al., this 
conference). The transverse width of the jet is of order of 0.5 — 0.7 mas, corre- 
sponding to a jet diameter of ~ 140 — 200 Schwarzschild radii (at r = 0.5 — 2 mas 
core separation). There are clear signs of edge-brightening and a 'hollow' or at 
least faint central jet body. We note that for a jet inclination > (30 — 45)°, a 
fast jet spine could be Doppler-deboosted and remain faint. 

In 2003, a component located ~ 0.2 mas east of a ~ 4 times brighter and 
more compact component is seen. At present it is unclear, if this eastern com- 
ponent should be identified with the jet base (which then must vary in com- 
pactness), or if the eastern component is part of a counter-jet. In any case, the 
structural differences between the 3 maps (and in the corresponding visibility 
data) indicate pronounced jet variability on and below the 1 year timescale. For 
the variations with an angular rate > 0.2mas/yr the corresponding jet speed 
would be /3app > 0.05. More frequent global 3mm- VLBI observations are re- 
quired to follow this motion. 
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6. VLBI at shorter wavelengths and Future Outlook 

After the successful detection of a number of radio source on continental and on 



transatlantic VLBI baselines at 129, 147 and 230 GHz (see lKrichbaum et al. 20071 . 
and references therein) , the technical feasibility of VLBI at wavelengths shorter 
than 3 mm has now been demonstrated. In the next step, and after upgrade of 
the existing VLBI data acquisition systems to higher data rates of ~ 2 — 4 Gbps 
and observing bandwidths of > 512 MHz, the baseline sensitivity will dramati- 
cally improve 0.5 Jy at 230 GHz for the existing antennas). With a relatively 
small VLBI array consisting of = 2 — 4 mm-telescopes, a systematic VLBI 
survey of the degree of compactness of a large sample of radio sources could be 
tackled. The further addition of more and more sensitive mm-capable telescopes 
(see Table [2|), will lead to true imaging capabilitj0. Such effort appears justi- 
fied, since the VLBI imaging of the Galactic Center source SgrA*, M87 and 
other nearby objects at 230 and 345 GHz will allow to map the direct vicinity 
of super massive black holes with a spatial resolution of < 5 — 6 Schwarzschild 
radii ! In this context, the participation of the large and sensitive millimeter 
interferometers (SMA, PdB, CARMA, ALMA) and large single dish telescopes 
(e.g the LMT) is essential, not only for reaching a high detection sensitivity, but 
also for providing a dense uv-coverage. Only with the latter an image fidelity 
can be reached, which in the end might allow to directly observe the expected 
image distortions and General Relativity effects in the vicinity of a black hole. 
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Name Country D Alt. Surf. 90 150 230 345 
[m] [m] [/xm] [GHz] [GHz] [GHz] [GHz] 



Surface < 50/xni 



SMA 


HI, USA 


8x6 


4100 


12 




+ 


+ 


HHT-SMTO 


AZ, USA 


10 


3100 


15 


V 


V 


+ 


APEX 


Chile 


12 


5000 


18 




+ 


+ 


ASTE 


Chile 


10 


4800 


20 


? 


? 


+ 


CSO 


HI, USA 


10 


4100 


25 


? 


+ 


+ 


JCMT 


HI, USA 


15 


4100 


25 


? 


V 


+ 



Surface 50 — lOO/im 



Pl.de Bure 


France 


6x15 


2550 


55 


V 


+ 


V 


+ 


Pico Veleta 


Spain 


30 


2900 


67 


V 


V 


V 


+ 


CARMA 


CA, USA 


a 


2200 


30-60 


+ 


? 


V 


+ 


Kitt Peak 


AZ, USA 


12 


2000 


75 


V 


V 


+ 





Surface 100 — 


300^lm 












Metsahovi 


Finland 


14 


sea 1. 


100 


V 


V 


Mopra 


Australia 


22 


900 


~120 


+ 




Onsala 


Sweden 


20 


sea 1. 


130 


V 


? 


VLBA 


USA 


25 


6 


150 


V 


? 


Taeduk 


Korea 


14 


200 


170 


+ 


? 


Nobeyama 


Japan 


45 


1300 


200 


+ 


+ 


Cambridge 


UK 


32 


20 


200 


? 




Haystack 


MA, USA 


37 


150 


250 


V 




Delingha 


China 


14 


3200 


250 


+ 




Simeiz 


Ukraine 


22 


sea 1. 


250 


? 





Surface 300 

GET 

Noto 

Effelsberg 



500/um 
VA, USA 
Italy 
Germany 



100 


840 


390^ 


+ 








32 


30 


400 


+ 








100 


300 


450 


V 









under construction 



ALMA 


Chile 


50x12 


5000 


25 


+ 


+ 


LMT 


Mexico 


50 


4600 


70 


+ 


+ 


Haystack*^ 


MA, USA 


37 


150 


100 


+ 


? 


Yebes 


Spain 


40 


850 


150 


+ 


? 


KVN 


Korea 


3x21 


200 


~150 


+ 


? 


SRT 


Italy 


64 


600 


~200 


+ 


? 



Notes: 

°: inhomogeneous array (6xl0m+9x6m) 

10 different antennas with elevations ranging from sea level to 4100 m 

present value, aim is ~ 210/um 
"^r after upgrade 

Table 2. Antennas for mm-VLBI, present status and future candidates. 
Col. 3 gives the antenna diameter, col. 4 the altitude, col. 5 the surface rms. 
Columns 6-9 summarize the suitability for mm-VLBI at a given observing 
band, with the following symbols: "V" : VLBI successfully done, receiver 
available; "+" : in principle possible, receiver available or planned; "?" : 
perhaps possible, but presently no receiver; "-" : not possible, surface not 
good enough. 



